I 


US  Army  Corps 
of  Engineers® 

Cold  Regions  Research  & 
Engineering  Laboratory 

Effect  of  Dissolved  NaCI  on 
Freezing  Curves  of  Kaolinite, 
Montmorillonite,  and  Sand  Pastes 

S.A.  Grant,  G.E.  Boitnott,  and  A.R.  Tice  January  1999 


pisrimmoy 

Appsov&d  ftss  pafoMe  tsS&aimi  f 
UaUxolt®dl  ' 


lyric  QUALITY  mSPECTED  1 


1 999030 1018 


Abstract:  We  developed  a  chemical-thermodynamic 
procedure  for  calculating  the  capillary  pressures  of  aque¬ 
ous  NaCI  solutions  in  a  porous  medium  at  tempera¬ 
tures  below  0°C  by  extending  the  treatment  by  Brun 
etal.  (1977).  Ice  in  the  porous  medium  was  assumed 
to  be  a  pure  phase  with  thermophysical  properties 
identical  to  bulk  hexagonal  ice.  The  thermophysical 
properties  (and  the  attendant  derivative  and  integral 
properties)  of  the  electrolyte  solutions  were  calculated 
with  the  Pitzer  model  as  parameterized  by  Archer 
(1992). 

Experiments  were  conducted  to  test  this  procedure. 
Pastes  of  kaolinite  clay,  montmorillonite,  and  quartz 


sand  were  prepared  by  washing  repeatedly  with  aque¬ 
ous  solutions  of  0.1  0.01  -,  and  0.001  -mol  kg-1  NaCI. 
The  molar  unfrozen  water  contents  of  these  pastes  were 
measured  by  pulsed  nuclear  magnetic  resonance  (NMR) 
in  the  temperature  range  -0.1 4°C  to  -66.6°C. 

The  relationships  between  ice-solution  capillary  pres¬ 
sures  and  specific  solution  volumes  for  frozen  pastes  of 
each  mineral  were  plotted  for  all  initial  solution  molali¬ 
ties.  While  some  systemic  errors  were  evident,  these  plots 
indicated  that  the  capillary  pressure-volume  relationships 
were  consistent  for  pastes  of  the  three  minerals  and,  as 
expected  from  theory,  unaffected  by  initial  equilibrating 
solution  molality. 
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Eq.  first 


Symbol 

Definition 

SI  unit  cited 

fl±NaCl(aq) 

mean-ionic  activity  of  NaCl(aq) 

dimensionless  59 

Vb 

activity  of  solute  B  (molality  concentration  scale) 

dimensionless  66 

ax,B 

activity  of  solute  B  (mole-fraction  concentration 

dimensionless  67 

scale) 

AC 

Debye-Huckel  coefficient  for  apparent 

J  K_1  mol-1 

46 

molar  constant-pressure  heat  capacity 

Alfm 

molar  area  of  solid/ gas  interface 

m2  mol-1 

3 

Als 

molar  area  of  liquid/ solid  interface 

m2moH 

3 

a 

fitted  coefficient  in  Maier-Kelly  equation 

J  K"1  mol"1 

33 

b 

fitted  coefficient  in  Maier-Kelly  equation 

J  K~2  mol"1 

33 

b 

a  constant  equal  to  1.2 

kg1/2  mol"1//2 

46 

Bp> 

^ p 

parameter  in  Speedy's  equation 

J  mol"1 

42 

c 

fitted  coefficient  in  Maier-Kelly  equation 

J  K  mol"1 

33 

^NaCl 

ion-interaction  parameter  of  Pitzer  model 

kg2  mol“2 

48 

r*(l) 

uNaCl 

ion-interaction  parameter  of  Pitzer  model 

kg2  mol”2 

48 

constant-pressure  heat  capacity 

JK-1 

42 

c 

m 

constant-pressure  molar  heat  capacity  of  the 

JK"1  mol"1 

44 

solution 

p0 

lp,h2o(1) 

constant-pressure  standard-state  heat 

J  K"1  mol"1 

44 

Cp(mr) 

capacity  of  liquid  water 
constant-pressure  heat  capacity  of  a 

J  K”1  mol"1 

46 

solution  with  a  molality  equal  to  mr 

Cq 

V 

parameter  in  Speedy's  equation 

J  mol"1 

42 

Im 

molality-based  ionic  strength 

mol  kg-1 

46 

lr 

molality-based  ionic  strength  of  the  reference 

mol  kg"1 

46 

solution 

m 

molal  concentration  of  the  solute 

mol  kg”1 

46 

«i 

initial  NaCl  molality  of  solutions  of  pastes 

mol  kgr1 

22 

mTf 

NaCl  molality  of  solutions  of  pastes  at  T=Tf 

mol  kg-1 

22 

mr 

reference  molality 

mol  kg-1 

46 

m±NaCl(aq) 

mean-ionic  molality  NaCl  in  an  aqueous 

mol  kg"1 

59 

solution 

m±NaCI(aq) 

standard-state  mean-ionic  molality  of  NaCl 

mol  kg"1 

59 

nT 

total  number  of  molecules  in  a  mixture  or 

mol 

36 

solution 

«A 

number  of  molecules  of  component  A  in 

mol 

36 

a  mixture  or  solution 

”H20(1) 

amount  of  water  in  the  solution 

mol 

44 

nNaCl(aq) 

amount  of  NaCl  in  the  solution 

mol 

44 

V 

pressure 

Pa 

1 

v 


pg 

pressure  of  gas  phase 

Pa 

i 

pi 

pressure  of  liquid  phase 

Pa 

i 

ps 

pressure  of  solid  phase 

Pa 

i 

Pc 

capillary  pressure 

Pa 

ii 

R 

gas  constant  (=  8.314  510  70) 

J  K-1  mol-1 

40 

Ss 

Dm 

molar  entropy  of  gas  phase 

J  K-1  moH 

1 

C*g 

molar  entropy  of  a  pure  water  vapor 

J  K”1  mol"1 

4 

Sm,H20 

molar  entropy  of  a  pure  water  liquid 

J  K-1  mol-1 

5 

Sm,H20 

molar  entropy  of  a  pure  water  ice 

J  K-1  mol"1 

6 

s1 

'-’m 

molar  entropy  of  liquid  phase 

J  K-1  mol"1 

1 

Ss 

Jm 

molar  entropy  of  solid  phase 

J  K-1  moH 

1 

T 

thermodynamic  temperature 

K 

1 

T0 

bulk  solvent  melting  point  for  a  given 
solution  composition 

K 

13 

Ts 

apparent  limiting  temperature  for 
supercooled  water  =  227.15  K 

K 

42 

V* 

v  m 

molar  volume  of  gas  phase 

m3 

1 

y1 

molar  volume  of  liquid  phase 

m3 

1 

Vs 

vm 

molar  volume  of  solid  phase 

m3 

1 

y*& 

Vh2o 

molar  volume  of  pure-water  vapor 

m3 

4 

T/*l 

vm/H20 

molar  volume  of  pure-water  liquid  phase 

m3 

5 

^m,H20 

molar  volume  of  pure-water  ice 

m3 

6 

*NaCl(aq) 

mole  fraction  of  B  in  the  a  phase 

dimensionless 

1 

zNa 

charge  number  of  Na+  cation 

dimensionless  46 

ZC1 

charge  number  of  Cl-  anion 

dimensionless  46 

a 

constant  in  Pitzer  model 

kg-l/2  mol-1/2  47 

a2 

constant  in  Pitzer  model  (as  revised  by  Archer) 

kg-1/2  mol”1/2  48 

P(NaCl 

ion-interaction  parameter  of  Pitzer  model 

kg  mol-1 

47 

P(NaCl 

ion-interaction  parameter  of  Pitzer  model 

kg  mol-1 

47 

AsHm,H20 

molar  enthalpy  of  melting  for  pure  water 

J  mol-1 

31 

y  Is 

liquid-solid  interfacial  tension 

Nm-1 

3 

ysg 

solid-gas  interfacial  tension 

Nm"1 

3 

Y±NaCl(aq) 

mean-ionic  activity  coefficient  for  NaCl  in  an 
aqueous  solution 

T  -Ts 

dimensionless  59 

e 

reduced  temperature  -  T 

1s 

dimensionless 

42 

HB 

chemical  potential  of  B  in  the  a  phase 

J  mol”1 

1 

M-m,B 

standard-state  chemical  potential  of  solute  B 
(molality  concentration  scale) 

J  mol”1 

60 

standard-state  chemical  potential  of  solute  B 
(mole-fraction  concentration  scale) 

J  mol-1 

61 

V 

=  vNa  +  VC1 

dimensionless 

46 

VI 


stoichiometric  number  of  Cl  in  NaCl  dimensionless  46 

stoichiometric  number  of  Na  in  NaCl  dimensionless  46 

initial  specific  liquid-water  content  kg  kg-1  22 

specific  liquid-water  content  at  T  =  Tf  kg  kg"1  22 

superscript  indicating  a  standard  substance  38 

superscript  indicating  a  pure  substance  4 

superscript  indicating  a  solution  at  infinite  dilution  60 


Effect  of  Dissolved  NaCl  on  Freezing  Curves  of 
Kaolinite,  Montmorillonite,  and  Sand  Pastes 

S.A.  GRANT,  G.E.  BOHNOTT,  AND  A.R.  TICE 


INTRODUCTION 

Chemical  thermodynamics  has  long  been  used  to  calculate  the  heaving  pres- 
sures  and  liquid-water  contents  of  frozen  ground  (Everett  1961,  Defay  and  Prigogine 
1966).  In  most  instances,  the  effects  of  solutes  in  the  liquid  phase  have  been  ne¬ 
glected  in  the  application  of  chemical-thermodynamic  theory  to  frozen  porous 
media:  the  liquid  phases  are  typically  assumed  to  be  pure.  At  equilibrium,  ice  ac¬ 
cepts  few  solutes  and  these  at  only  vanishingly  small  mole  fractions.  Most  solutes 
in  frozen  ground  are  in  the  liquid  phase  (Hobbs  1974).  At  temperatures  apprecia¬ 
bly  below  0°C,  where  most  of  the  water  in  a  water-electrolyte  system  is  in  the  ice 
phase,  liquid  solutions  tend  to  be  highly  concentrated.  In  these  environments,  the 
effects  of  solutes  on  the  freezing  behavior  of  ground  can  be  considerable. 

In  recent  years,  theoretical  and  experimental  advances  have  provided  the  tools 
for  estimating  the  thermophysical  properties  of  aqueous  electrolyte  solutions  at 
subzero  temperatures.  The  major  theoretical  advance  has  been  the  development 
and  acceptance  of  the  Pitzer  model  of  electrolyte  solutions  (Pitzer  1991).  The  Pitzer 
model  allows  for  the  calculation  of  the  important  physiochemical  properties  of 
simple  and  complex  aqueous  electrolyte  solutions  over  a  wide  range  of  tempera¬ 
tures,  pressures,  and  compositions.  Most  importantly  for  this  discussion,  the  Pitzer 
model  allows  the  accurate  estimation  of  electrolyte-solution  properties  at  subzero 
temperatures  (Spencer  et  al.  1990,  Archer  1992,  Marion  and  Grant  1994).  Since  many 
of  the  thermophysical  properties  modelled  by  the  Pitzer  model  are  relative  to  the 
same  properties  of  the  solvent,  a  critical  component  when  modelling  the 
thermophysical  properties  of  electrolyte  solutions  is  accurate  measurements  of  pure- 
solvent  properties.  Historically,  this  has  presented  problems  for  water  at  tempera¬ 
tures  below  its  equilibrium  freezing  point.  Novel  experimental  procedures  have 
allowed  the  heat  capacity  and  density  of  supercooled  liquid  water  to  be  measured 
accurately  at  temperatures  as  low  as  -35°C  (Angell  et  al.  1973,  Hare  and  Sorensen 
1987).  Speedy  (1987)  has  fitted  much  of  the  available  data  on  supercooled  water  to 
a  family  of  equations  that  calculate  physical-chemical  properties  of  liquid  water 
under  atmospheric  pressure  in  the  temperature  range  0  to  -45°C. 

To  our  knowledge,  no  one  has  applied  these  new  capabilities  to  the  long-stand¬ 
ing  problem  of  calculating  the  effects  of  solutes  on  the  freezing  curves  of  frozen 
ground.  The  objectives  of  the  research  presented  here  were: 

1.  To  extend  the  conventional  (i.e.,  pure-phase)  chemical-thermodynamic  theory 
of  frozen  ground  to  the  liquid  phases  of  aqueous  electrolyte  solutions 

2.  To  calculate  the  needed  chemical-thermodynamic  data  from  physical-chemi¬ 
cal  models  of  ice,  water,  and  electrolyte  solutions 

3.  To  make  experimental  measurements  of  liquid-water  contents  of  sand  and 
clay-mineral  pastes  with  which  to  test  this  extension  of  the  model. 

In  this  report,  we  initially  discuss  the  chemical-thermodynamic  methods  by 
which  capillary  pressures  of  either  pure  or  solution  liquid  phases  can  be  calculated 
for  frozen  porous  media.  The  thermophysical  data  needed  to  make  these  calcula¬ 
tions  are  presented.  We  then  describe  and  evaluate  the  data  from  a  set  of  experi¬ 
ments  we  conducted  to  test  this  approach  to  calculating  capillary  pressures. 


THEORY 


Here  we  follow  closely  the  lucid,  complete  presentation  of  the  chemical-ther¬ 
modynamic  theory  of  frozen  porous  media  presented  by  Brun  et  al.  (1977).  Con¬ 
sider  three  phases,  solid  (s),  liquid  (1),  and  gas  (g),  which  are  in  thermal,  but  not 
necessarily  hydrostatic,  equilibrium  [i.e.,  T8  =  Tl  =  TS  =  T,  but  ps  *  pl  *  p,  where  T  is 
temperature  (K)  and  p  is  pressure  (Pa)].  Their  development  begins  with  two  sets  of 
classic  equalities: 

1.  The  Gibbs-Duhem  equations,  and 

2.  The  Kelvin  equations. 

Gibbs-Duhem  equations 

For  each  phase,  the  Gibbs-Duhem  relation  holds: 

o  =  SmdT  -  Vgdpg  + 1  xfdn! 

B 

0  =  SimdT-Vlmdpl+'L*WB  (1) 

0  =  S^dT  -  v^dps  + 1  *|d|i| 

B 

where  S^,  S^,  and  are  the  molar  entropies  (J  K-1  mol-1)  of  the  gas,  liquid,  and 
solid  phases,  respectively;  V^,  V^,  and  V£,  their  molar  volumes  (m3  mol-1);  x\, 
jtg,  and  x%,  the  mole  fraction  of  component  B  in  the  three  phases;  and  jli|,  jig,  and 
p|,  the  chemical  potentials  (J  mol-1)  of  component  B  in  the  three  phases. 

Kelvin  equations 

The  Kelvin  equations  give  the  pressure  gradients  across  the  gas/ solid  and  liq¬ 
uid/solid  interfaces  in  terms  of  the  differential  geometry  of  these  interfaces: 


d As% 

pS_ps=yssri^EL 

dVg 


(2) 


pl_ps=Yb^k 


dV, 


(3) 


where  and  A*fm  (m2  mol-1)  are  the  molar  areas  of  the  liquid/solid  and  solid/ 
gas  interfaces,  respectively,  and  yls  and  y  SS  (N  m-1)  are  the  corresponding  interfa¬ 
cial  tensions. 

Capillary  pressures  of  liquid  phases  can  be  estimated  by  applying  these  equali¬ 
ties.  Below  we  present  a  summary  of  the  development:  first  for  systems  with  pure- 
water  liquid  phases,  then  for  aqueous-electrolyte-solution  liquid  phases. 


Systems  with  pure  liquid-water  phases 

If  the  solid,  liquid,  and  vapor  phases  are  pure,  then  eq  1  can  be  rewritten: 


0  =  Sm,H2OdT  “  Vm,K2OdPS  +  d^H20 

(4) 

0  =  Sm,H2Odr  “  Vm,H2OdP1  +  dM-H20 

(5) 

0  =  Sm,H2OdT  -  ^irf,H2OdPS  +  d^H20 

(6) 

where  the  superscript  *  indicates  a  pure  phase.  By  taking  the  partial  derivatives  of 
eq  2  and  3,  dps  and  dp1  can  be  related  to  dps  by: 


2 


(7) 


(  cL4sg 

dpS=dps+d  YSS— J* 


and 


V 


dp1  =  dps+d  yk-— ^ 
{  dv* 


(8) 


By  subtracting  eq  4  from  eq  6  and  eq  6  from  eq  5,  then  subtracting  the  latter  differ¬ 
ence  from  the  former,  the  following  relation  is  obtained  after  rearrangement: 


Sm,H2Q  -Sm,H2Q 

v*s  t/S 
vm,H20  “  Vm,H20 


Sm,H2Q 

T/*l 

vm,H20 


c*s 

T/+S 

“  '/m/H20 


\ 

dT 

y 


y*S 

V.HzO 


VmgH20  -  Vm,H20 


7  cMSg 
sg  gjVn 

dVS  , 


(9) 


T/*l 

*m,H20 


Vm,H20  “  ^m,H20 


1 


V 


dV, 


m  ) 


By  noting  that  V^h2o  »  and  by  assuming  that  ysg  — =  0,  the  follow¬ 
ing  relation  is  obtained:  dlm 


(Sm,H20  “  Sm,H2o)dT  =  ^HjOd 


dym  , 


By  defining  capillary  pressure  (Pa),  pc,  as 


(10) 


Pc=P1-PS' 


(11) 


it  can  be  seen  from  eq  8  and  11  that 


dPc=d^ 


cL4f 


dV* 


A 


J 


(12) 


Therefore,  the  capillary  pressure  of  liquid  water  in  a  frozen  porous  medium  can  be 
calculated  by 


Prcj  Ti  Sm,H20-Sm,H20 

jdpc=pc=  j - Pi - — 

0  To  '/m,H20 


dT 


(13) 


where  T0  (K)  is  the  melting  point  of  the  bulk  liquid  phase. 

To  evaluate  the  righthand  side  of  eq  13,  the  following  properties  must  be  mea¬ 
sured  or  estimated  for  all  temperatures  of  interest: 

1.  The  melting  point  of  the  liquid  phase  in  bulk 

2.  Entropy  of  the  liquid  phase 

3.  Entropy  of  the  ice  phase 

4.  Molar  volume  of  the  liquid  phase. 
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Systems  with  aqueous-electrolyte-solution  liquid  phases 

The  treatment  presented  in  the  previous  section  can  be  extended  directly  to  elec- 
trolyte-H^O  systems.  For  example,  eq  1  for  a  NaCl-H20  system  can  be  written 


0  =  CgH2OdT-VmgH2Od^+d^H2O 

(14) 

o  =  SjndT  -  v}ndp  +  *H2od(iH2o  +  XNaCldflNaCl 

(15) 

0  =  Sm,H2odT  “  VmS,H2OdPS  +dltH20 

(16) 

For  a  simple  NaCl  aqueous  solution,  x{sjaci  =  1  -  *h2o  following  relation  is  de¬ 

rived  by  steps  similar  to  those  that  led  to  eq  9: 


^  c*s  c*S 

bm,H2Q  ~  ^m,H2Q 

^ ^m,H20  ”  ^m,H20 


\ 

■  dT 
) 


^mgHzO  d 

cL4sg  "l 

sg 

V1 

m  d 

d/4ls  ^ 
..Is  “As.m 

ymgH20_VmS,H20 

y 

iVg  ) 

“  ^H20^mS,H20 

did 

^  m  y 

*NaCl 


^m“4i2O^mS,H20 


^M'NaCl 


(17) 


d  A* 


As  before,  by  noting  that  ^gH2o  »  ym,H20  and  by  assuming  that  ysg  S/™  =  0, 
the  following  relation  is  obtained:  m 

(  cL4ls  1 

(sJn  -  ^2oS^,h2o)  dT  =  V^d  Yls  — Sf  -  -  -’•NaCld^NaCl  (18) 


which  yields  on  rearrangement 


vid 


(  H/*ls  ^  ( 

Is 

dVl 


Sm  “  *H,0Sm,H,0  +  xNaCl 


d^NaCl 

dT 


dT  . 


(19) 


The  ice-solution  capillary  pressure  can  be  calculated  as  a  function  of  temperature 
with 


rl  ..1  r*S  .  ..l  dpNaCl 

m  -  xH20i,m,H20  +  *NaCl  — 

Pc  =  \ - - dT- 


(20) 


rm 


To  calculate  ice-solution  capillary  pressure,  the  quantities  in  the  integrand  must  be 
calculated.  These  quantities  are 

1.  The  mole  fractions  of  the  solute  and  solvent 

2.  The  melting  point  of  the  solution  phase  in  bulk 

3.  The  molar  entropy  of  ice 

4.  The  molar  entropy  of  the  liquid  solution 

5.  The  molar  volume  of  the  liquid  solution 

6.  The  temperature  derivative  of  the  solute  chemical  potential. 
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Relationship  of  capillary  pressures  to  pore  radii  and  liquid-water  contents 

Defay  and  Prigogine  (1966)  have  derived  a  form  of  the  Laplace  equation  for 
freezing  liquids  in  saturated  porous  media: 


Pc=~ 


(21) 


where  rls  is  the  radius  of  a  concave  liquid /solid  interface  in  a  pore.  Brun  et  al. 
(1977)  found  that  rls  is  equal  to  the  pore  radius  less  the  few  molecular  diameters 
corresponding  to  the  liquid  that  wets  the  pore  solid  surfaces,  which  remains 
unfrozen  to  temperatures  well  below  the  bulk  melting  point.  Experimentally  mea¬ 
sured  freezing  curves  of  porous  solids  are  consistent  with  a  general  rela¬ 
tionship  between  liquid/solid  capillary  pressure  and  the  specific  volume 
of  unfrozen  liquid.  For  a  given  capillary  pressure  corresponding  to  a  particular 
pore  radius,  it  is  assumed  generally  that  the  contents  of  all  narrower  pores  are 
unfrozen  and  the  contents  of  all  wider  pores  are  solidified.  This  capillary 
pressure-unfrozen  liquid  relationship  is  physical  and  should  be  unaffected  directly 
by  the  chemical  composition  of  the  phases.  For  a  given  porous  material,  there¬ 
fore,  it  is  expected  that  this  relationship  will  be  general  and  unaffected  by 
phase  composition.  We  tested  this  expectation  with  experiments  described  in  the 
following  section. 


MATERIALS  AND  METHODS 
Sample  preparation 

Kaolinite  clay 

Approximately  40  g  of  clay  were  washed  five  times  with  200  mL  of  1-mol  kgr1 
NaCl  aqueous  solution  to  saturate  the  cation  exchange  complex  with  Na+.  In  each 
washing,  the  clay-solution  suspension  was  stirred  for  about  5-15  minutes  with  a 
vortex  stirrer.  The  solids  were  separated  with  a  centrifuge. 

Following  the  last  washing  with  1-mol  kg”1  NaCl  solution,  the  suspension  was 
separated  into  three  aliquots  and  individually  washed  with  0.1-,  0.01-,  and  0.001- 
mol  kg"1  NaCl  solutions,  respectively.  These  washings  were  followed  by  centrifu¬ 
gation;  this  process  was  repeated  four  times.  The  individual  samples  were  then 
sealed  in  test  tubes  with  rubber  stoppers  to  prevent  moisture  changes. 

Montmorillonite  clay 

Sodium-saturated  montmorillonite  that  had  been  prepared  previously  for 
other  experiments  and  then  stored  was  selected  for  analysis.  Three  4-g 
samples  of  this  clay  were  placed  in  beakers  to  which  200-mL  solutions  of 
either  0.1-,  0.01-,  or  0.001-mol  kg”1  NaCl  were  added.  The  clays  were  then  sepa¬ 
rated  from  the  solutions  and  resuspended  five  times  with  further  washings  of  200- 
mL  solutions  of  0.1-,  0.01-,  and  0.001-mol  kg”1  NaCl.  The  0.1-mol  kg”1  suspension 
was  centrifuged  at  high  speed  for  approximately  6  hours.  The  liquids  in  the  0.01- 
and  0.001-mol  kg-1  suspensions  were  separated  from  the  clay  by  pressure-mem¬ 
brane  extraction.  Following  the  above  treatments,  the  samples  were  placed  in  sealed 
test  tubes. 

Quartz  sand 

A  fine-textured  sand  was  divided  into  three  50-g  subsamples  that  were  washed 
four  times  with  either  0.1-,  0.01-,  or  0.001-mol  kg-1  NaCl  solution.  Each  saturated 
sand  sample  was  sealed  in  test  tubes. 
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Experimental  procedure 

Once  the  mineral  pastes  were  placed  in  their  respective  test  tubes,  at  least  two 
days  were  allowed  for  moisture  equilibration  before  the  unfrozen  water  content 
was  determined.  Following  this  period,  the  sealed  tubes  containing  minerals  and 
solutions  were  placed  in  a  bath  preset  to  an  initial  test  temperature  of  about  20°C.  If 
the  analysis  temperatures  were  within  the  range  of  20°C  to  -20°C,  the  bath  heat- 
transfer  fluid  was  an  ethylene  glycol-water  mixture.  If  the  analysis  required  tem¬ 
peratures  below  -20°C,  a  second  bath  that  used  methyl  alcohol  as  a  heat-transfer 
fluid  and  was  capable  of  achieving  temperatures  of  -70°C  was  utilized. 

Following  a  temperature  equilibration  time  of  approximately  1  hour,  the  sample 
tubes  were  sequentially  removed  from  the  bath,  wiped  dry,  and  inserted  into  the 
nuclear  magnetic  resonance  (NMR)  probe  for  analysis.  (The  mean  background  sig¬ 
nal  intensity  due  to  the  test  tube  was  measured  by  placing  an  empty  test  tube 
in  the  NMR  probe.)  The  peak  intensity  for  each  temperature  was  recorded.  The 
NMR  probe  used  for  these  analyses  was  a  Praxis  model  PR-103  analyzer.  It  was  oper¬ 
ated  in  the  90°  mode  with  a  0.2-s  clock  and  at  a  fast  scan  speed.  The  first  pulse  ampli¬ 
tude  in  the  90°  mode  was  measured  for  each  sample,  starting  with  the  first  test  tem¬ 
perature  of  about  20°C.  After  about  4  s  (the  time  required  to  measure  the  NMR  signal 
amplitude)  the  samples  were  reinserted  into  the  bath.  When  all  the  samples  were 
analyzed,  the  bath  temperature  was  lowered  by  about  3°C  and  allowed  to  re-equili- 
brate  at  this  new  test  temperature.  This  process  was  repeated  until  a  temperature  of 
about  0°C  was  reached.  (All  of  the  above-0°C  NMR  measurements  were  used  to  de¬ 
termine  the  paramagnetic  effect  or  the  increase  in  signal  intensity  with  decreasing 
temperature.  This  relationship  is  discussed  by  Tice  et  al.  [1981, 1982].) 

Once  measurements  above  0°C  were  completed,  the  constant-temperature  bath 
was  cooled  to  about  -5°C  to  initiate  ice  formation  in  the  samples.  The  constant- 
temperature  bath  was  then  raised  to  -0.5°C,  and  the  samples  were  allowed  to  equili¬ 
brate  at  this  temperature  for  several  hours.  After  equilibration,  the  peak  intensity  and 
temperature  were  recorded.  The  temperature  of  the  constant-temperature  bath  was 
then  lowered  in  small  increments.  The  signal  intensity  and  temperature  were  recorded 
after  the  sample  was  allowed  to  equilibrate  at  each  temperature  increment.  The  samples 
were  cooled  until  the  measured  signal  intensity  matched  that  for  an  empty  test  tube — 
indicating  that  no  liquid  water  remained  in  the  sample  (Tice  1982). 

Following  the  last  NMR  cooling  measurements,  which  in  some  instances  were 
taken  at  temperatures  as  cold  as  -77°C,  the  warming  run  was  begun.  An  analysis 
interval  of  about  4°C  was  selected  at  the  lower  temperatures,  at  which  small  changes 
in  temperature  have  little  effect  on  signal  intensity.  At  temperatures  near  melting, 
intervals  between  equilibration  temperatures  were  reduced  because  liquid-water 
contents  were  increasing  more  with  temperature.  At  around  -1°C,  the  interval  by 
which  the  temperature  was  changed  was  reduced  to  about  -0.1°C.  This  approach 
provided  good  coverage  within  the  region  where  the  unfrozen  water-temperature 
relationship  is  so  critical. 

Unfrozen  water  contents  were  calculated  by  first  regressing  the  above-0°C  NMR 
readings  minus  the  background  against  their  respective  temperatures  for  each 
sample  and  extending  the  resultant  line  to  the  lowest  temperature  where  experi¬ 
mental  data  were  obtained.  Projected  first-pulse  amplitudes  were  calculated  for 
each  experimental  temperature.  The  gravimetric  water  content  for  each  sample  had 
previously  been  determined  by  oven-drying  at  110°C  for  two  days.  A  ratio  between 
the  sample  water  content  and  the  projected  first-pulse  amplitude  was  developed; 
unfrozen  water  contents  were  calculated  by  multiplying  these  first-pulse  ampli¬ 
tudes  by  their  respective  ratios  to  obtain  a  value  for  each  temperature  (Tice  et  al. 
1981, 1982). 
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EXPERIMENTAL  RESULTS 


The  mean  initial  specific  water  contents  and  equilibrating-solution  molalities 
for  the  pastes  are  presented  in  Table  1.  The  specific  liquid  water  contents  of  the 
pastes  as  measured  by  pulsed  NMR  are  presented  in  Tables  2  through  7,  and  their 
liquid  water  contents  are  presented  graphically  in  Figures  1  through  6. 

For  a  given  initial  solution  molality,  the  freezing  curves  of  all  pastes  displayed 
the  roughly  exponential  decline  in  liquid-water  content  with  decreasing  tempera¬ 
ture  that  is  characteristic  of  freezing  curves.  The  freezing  curves  of  the  sand  pastes 
were  much  less  smooth  than  those  for  the  kaolinite  or  montmorillonite  pastes  be¬ 
cause  of  the  much  lower  intensity  of  the  peak  signal  from  the  sand  samples.  There 
was  no  pronounced  hysteresis  between  the  cooling  and  heating  curves  for  any  of 
the  samples.  For  each  of  the  three  minerals,  the  freezing  curves  of  the  pastes  that 
had  been  equilibrated  with  0.001-mol  kg-1  NaCl  solutions  were  roughly  coincident 
with  those  that  were  equilibrated  with  0.01-mol  kg^1  NaCl  solutions.  For  each  min¬ 
eral,  the  freezing  curves  of  the  pastes  initially  equilibrated  with  0.1-mol  kg-1  NaCl 
solutions  showed  more  liquid  water  at  a  given  temperature  than  those  equilibrated 
with  the  lower-molality  solutions.  This  shifting  in  the  curves  was  likely  the  result 
of  the  depression  of  the  freezing  points  of  the  higher-molality  solution.  On  a  rela¬ 
tive  basis,  this  shift  was  most  pronounced  for  the  sand  samples  and  least  pronounced 
for  the  kaolinite  samples. 

COMPARISON  WITH  THEORY 

In  this  section  we  first  detail  how  the  capillary  pressures  and  specific  volumes  of 
the  pore  water  solutions  were  calculated.  These  calculated  values  are  then  pre¬ 
sented  graphically  and  analyzed. 

Calculation  of  the  necessary 
thermophysical  properties 

As  stated  earlier,  the  following  thermodynamic  quantities  are  needed  to  calcu¬ 
late  ice-solution  capillary  pressures  in  frozen  porous  media: 

1.  The  mole  fractions  of  the  solute  and 
solvent 

2.  The  freezing  point  of  the  solution 
phase  in  bulk 

3.  The  molar  entropy  of  ice 

4.  The  molar  entropy  of  the  liquid  solu¬ 
tion 

5.  The  molar  volume  of  the  liquid  solu¬ 
tion 

6.  The  temperature  derivative  of  the  sol¬ 
ute  chemical  potential. 

Mole  fractions  of  the  solute  and  solvent 

A  bulk  H20-NaCl  system  under  atmo¬ 
spheric  pressure  has  a  precisely  defined  eu¬ 
tectic  point  at  252  K  (-21.15°C)  and  5.14  mol 
kg-1.  Below  the  eutectic  temperature  no  liq¬ 
uid  phase  remains,  giving  way  to  a  mixed 
solid  phase  composed  of  hexagonal  ice  [rep- 


Table  1.  Equilibrating  solution  mo¬ 
lalities  and  total  specific  water 
masses  for  pastes  frozen  then 
thawed  in  this  study. 


Mineral 

Initial 

NaCl 
molality 
(mol  kg-1) 

Initial 

specific  water 
content 
(kg  kg'1) 

Kaolinite 

0.100  471 

2.3425 

0.010  031  4 

2.6703 

0.001  002  98 

2.1865 

Montmoril¬ 

0.100  471 

35.8610 

lonite 

0.010  031  4 

31.5064 

0.001  002  98 

12.5529 

Sand 

0.100  471 

0.2165 

0.010  031  4 

0.2146 

0.001  002  98 

0.2225 
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Table  2.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  kaolin- 
ite  pastes  cooled  from  0°C  to  -66.6°C. 

Specific  liquid  water  content  (kg  kg'1) 
Initial  NaCl  solution  molality 


Temperature 

VC) 

0.1 

mol  kg1 

0.01 
mol  kg1 

0.001 
mol  kg~J 

-0.25 

0.9873 

0.7821 

-0.34 

0.8281 

0.6542 

-0.43 

0.6231 

0.5024 

-0.52 

1.6059 

0.5011 

0.4332 

-0.61 

1.5677 

0.4769 

0.4255 

-0.76 

1.3924 

0.4361 

0.3982 

-0.81 

1.2570 

0.4267 

0.3817 

-0.89 

1.1216 

0.4026 

0.3755 

-1.10 

0.8735 

0.3838 

0.3721 

-1.20 

0.8146 

0.3780 

0.3480 

-1.29 

0.7493 

0.3650 

0.3419 

-1.45 

0.6774 

0.3482 

0.3371 

-1.54 

0.6333 

0.3498 

0.3325 

-1.76 

0.5856 

0.3439 

0.3083 

-1.96 

0.5284 

0.3124 

0.2737 

-2.15 

0.4875 

0.2719 

0.2437 

-2.43 

0.4562 

0.2478 

0.2211 

-2.69 

0.4233 

0.2111 

0.1986 

-3.02 

0.4017 

0.1998 

0.1790 

-3.21 

0.3868 

0.1723 

0.1537 

-3.75 

0.3553 

0.1501 

0.1371 

-4.40 

0.3046 

0.1135 

0.1072 

-5.70 

0.2328 

0.0841 

0.0803 

-7.20 

0.1917 

0.0656 

0.0579 

-9.00 

0.1416 

0.0474 

0.0431 

-11.00 

0.1063 

0.0347 

0.0327 

-14.10 

0.0773 

0.0273 

0.0238 

-17.00 

0.0598 

0.0219 

0.0207 

-20.00 

0.0515 

0.0165 

0.0163 

-23.00 

0.0421 

0.0130 

0.0107 

-24.70 

0.0432 

0.0145 

0.0120 

-26.70 

0.0356 

0.0080 

0.0092 

-28.80 

0.0324 

0.0079 

0.0117 

-30.60 

0.0308 

0.0063 

0.0104 

-33.80 

0.0276 

0.0077 

0.0089 

-39.00 

0.0188 

0.0075 

0.0062 

-45.20 

0.0052 

0.0044 

0.0024 

-49.80 

0.0038 

0.0048 

-55.30 

0.0025 

0.0047 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1, 0.01,  and  0.001  mol  kg*1  NaCl. 
Data  are  not  presented  for  temperatures  at  which 
no  liquid  water  was  detected. 


Table  3.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  kaolin- 
ite  pastes  warmed  from  -66.6°C  to  0°C. 

Specific  liquid  water  content  (kg  kg'1) 
Initial  NaCl  solution  molality 


Temperature 

(V) 

0.1 

mol  kg"1 

0.01 
mol  kg"1 

0.001 
mol  kg-1 

-61.20 

0.0024 

0.0023 

-55.60 

0.0013 

0.0028 

0.0012 

-50.10 

0.0038 

0.0029 

0.0036 

-44.80 

0.0052 

0.0073 

0.0061 

-39.30 

0.0081 

0.0090 

0.0087 

-34.60 

0.0096 

0.0092 

0.0114 

-30.60 

0.0182 

0.0078 

0.0117 

-28.20 

0.0212 

0.0127 

0.0157 

-26.40 

0.0314 

0.0144 

0.0132 

-24.50 

0.0389 

0.0129 

0.0160 

-22.40 

0.0436 

0.0131 

0.0121 

-19.20 

0.0532 

0.0183 

0.0177 

-17.30 

0.0597 

0.0235 

0.0193 

-13.60 

0.0790 

0.0257 

0.0224 

-10.40 

0.1051 

0.0383 

0.0342 

-7.90 

0.1534 

0.0530 

0.0519 

-6.00 

0.1992 

0.0732 

0.0627 

-4.20 

0.2874 

0.1064 

0.0926 

-3.26 

0.3578 

0.1578 

0.1389 

-2.80 

0.3877 

0.1782 

0.1644 

-2.18 

0.4649 

0.2391 

0.2006 

-1.75 

0.5226 

0.2781 

0.2517 

-1.63 

0.5537 

0.3149 

0.2787 

-1.50 

0.5767 

0.3316 

0.2863 

-1.33 

0.6113 

0.3558 

0.3209 

-1.20 

0.6377 

0.3670 

0.3316 

-1.13 

0.6672 

0.3764 

0.3511 

-1.01 

0.6855 

0.3895 

0.3618 

-0.97 

0.7003 

0.3657 

0.3499 

-0.86 

0.7332 

0.3806 

0.3501 

-0.80 

0.7708 

0.4010 

0.3637 

-0.65 

0.8040 

0.4142 

0.3715 

-0.58 

0.8353 

0.4384 

0.3896 

-0.50 

0.8780 

0.5012 

0.4123 

-0.34 

0.9293 

0.5939 

0.4876 

-0.27 

0.9851 

0.6532 

0.5448 

-0.21 

0.7365 

0.6186 

-0.19 

0.7642 

0.6502 

-0.14 

0.8217 

0.6999 

-0.08 

0.8737 

0.7467 

-0.02 

0.9498 

0.8056 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1,  0.01,  and  0.001  mol  kg-1  NaCl. 
Data  are  not  presented  for  temperatures  at 
which  no  liquid  water  was  detected. 
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Table  4.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  mont- 
morillonite  pastes  cooled  from  0°C  to 
-66.6°C 

Specific  liquid  water  content  (kg  kg-1) 
Initial  NaCl  solution  molality 


Temperature 

CO 

0.1 

mol  kg'1 

0.01 
mol  kg-1 

0.001 
mol  kg'1 

-0.14 

2.5987 

1.8358 

-0.25 

1.9654 

1.5490 

-0.34 

1.6137 

1.1300 

-0.43 

1.2097 

0.8856 

-0.52 

1.1391 

0.6413 

-0.61 

0.8583 

0.5854 

-0.76 

0.7351 

0.5433 

-0.81 

0.7174 

0.5292 

-0.89 

0.6647 

0.4316 

-1.10 

0.6465 

0.4313 

-1.20 

0.5938 

0.3964 

-1.29 

0.7332 

0.3823 

-1.45 

0.6453 

0.3752 

-1.54 

0.6974 

0.3889 

-1.76 

0.6617 

0.3748 

-1.96 

0.6089 

0.3606 

-2.15 

0.5735 

0.3188 

-2.43 

16.0658 

0.5207 

0.3254 

-2.69 

5.3741 

0.5200 

0.3251 

-3.02 

4.6635 

0.5018 

0.3040 

-3.21 

4.2586 

0.4840 

0.2969 

-3.75 

3.8469 

0.4655 

0.2757 

-4.40 

3.3165 

0.3952 

0.2682 

-5.70 

2.6221 

0.3244 

0.2259 

-7.20 

2.0724 

0.2712 

0.1906 

-9.00 

1.5292 

0.1848 

0.1489 

-11.00 

1.1680 

0.1497 

0.1277 

-14.10 

0.8850 

0.1147 

0.1130 

-17.00 

0.7046 

0.1132 

0.1054 

-20.00 

0.5659 

0.0797 

0.1043 

-23.00 

0.4494 

0.0472 

0.1033 

-24.70 

0.4277 

0.0624 

0.0834 

-26.70 

0.4058 

0.0773 

0.0765 

-28.80 

0.0613 

0.0760 

-30.60 

0.0693 

-33.80 

0.0623 

-39.00 

0.0490 

-45.20 

0.0541 

-49.80 

0.0297 

-55.30 

0.0175 

-60.70 

0.0230 

-61.20 

-66.60 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1, 0.01,  and  0.001  mol  kg-1 
NaCl.  Data  are  not  presented  for  temperatures 
at  which  no  liquid  water  was  detected. 


Table  5.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  mont- 
morillonite  pastes  warmed  from 
-66.6°C  to  0°C. 

Specific  liquid  water  content  (kg  kg-1) 
Initial  NaCl  solution  molality 


Temperature 

ra 

-  0.1 
mol  kg'1 

0.01 
mol  kg'1 

0.001 
mol  kg1 

-61.20 

0.0000 

0.0000 

0.0000 

-55.60 

0.0000 

0.0000 

0.0350 

-50.10 

0.0000 

0.0000 

0.0415 

-44.80 

0.0000 

0.0000 

0.0361 

-39.30 

0.0332 

0.0293 

0.0613 

-34.60 

0.0509 

0.0448 

0.0559 

-30.60 

0.0519 

0.0760 

0.0756 

-28.20 

0.1051 

0.0768 

0.0825 

-26.40 

0.0883 

0.0619 

0.0830 

-24.50 

0.1249 

0.0937 

0.0963 

-22.40 

0.1262 

0.0788 

0.0970 

-19.20 

0.5682 

0.0800 

0.1046 

-17.30 

0.6109 

0.0807 

0.1052 

-13.60 

0.8118 

0.1150 

0.1132 

-10.40 

1.1910 

0.1668 

0.1347 

-7.90 

1.7528 

0.2533 

0.1494 

-6.00 

2.1847 

0.2898 

0.1573 

—1.20 

3.0220 

0.4300 

0.1858 

-3.26 

3.8376 

0.5185 

0.2761 

-2.80 

4.1884 

0.6410 

0.2904 

-2.18 

5.1082 

0.6603 

0.3188 

-1.75 

6.1489 

0.7314 

0.3609 

-1.63 

7.0005 

0.7145 

0.3819 

-1.50 

7.4095 

0.7498 

0.3959 

-1.33 

8.6090 

0.8377 

0.3753 

-1.20 

10.1929 

0.9081 

0.3825 

-1.13 

11.4315 

0.9434 

0.3895 

-1.01 

12.4109 

0.9090 

0.4384 

-0.97 

12.9605 

1.0665 

0.4385 

-0.86 

14.0224 

1.1021 

0.4804 

-0.80 

15.2435 

0.9450 

0.5084 

-0.65 

16.8189 

1.0333 

0.5783 

-0.58 

18.4088 

1.2439 

0.5854 

-0.50 

20.1838 

1.3320 

0.6762 

-0.34 

22.5804 

1.7716 

0.9068 

-0.27 

25.3344 

2.0179 

1.1093 

-0.21 

0.0000 

1.9659 

1.3189 

-0.19 

0.0000 

2.1240 

1.4307 

-0.14 

0.0000 

2.7216 

1.7939 

-0.08 

0.0000 

3.4426 

2.1993 

-0.02 

0.0000 

4.4276 

2.7235 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1,  0.01,  and  0.001  mol  kg-1 
NaCl.  Data  are  not  presented  for  temperatures 
at  which  no  liquid  water  was  detected. 
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Table  6.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  sand 
pastes  cooled  from  0°C  to  -66.6°C. 

Specific  liquid  water  content  (kg  kg'1) 
Initial  NaCl  solution  molality 
Temperature  0.1  0.01  0.001 

( <C)  mol  kg'1  mol  kg'1  mol  kg'1 


-0.14 

0.0145  0.0119 

-0.25 

0.0106  0.0069 

-0.34 

0.0096  0.0040 

-0.43 

0.0096 

-0.52 

0.0086 

-0.61 

0.0067 

-0.76 

0.0067 

-0.81 

0.0038 

-0.89 

0.0048 

-1.10 

0.0076 

-1.20 

0.0067 

-1.29 

0.0057 

-1.45 

0.0047 

-1.54 

0.0085 

-1.76 

0.0085 

-1.96 

0.0066 

-2.15 

0.0481 

0.0066 

-2.43 

0.0461 

0.0066 

-2.69 

0.0431 

0.0056 

-3.02 

0.0392 

0.0047 

-3.21 

0.0410 

0.0056 

-3.75 

0.0380 

0.0046 

^.40 

0.0340 

0.0046 

-5.70 

0.0290 

0.0045 

-7.20 

0.0214 

-9.00 

0.0149 

-11.00 

0.0111 

-14.10 

0.0075 

-17.00 

0.0056 

-20.00 

0.0039 

-23.00 

0.0023 

-24.70 

0.0023 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1, 0.01,  and  0.001  mol  kg-* 
NaCl.  Data  are  not  presented  for  temperatures 
at  which  no  liquid  water  was  detected. 

Table  7.  Unfrozen-water  contents,  as 
measured  by  pulsed  NMR,  of  sand 
pastes  warmed  from  -66.6°C  to  0°C. 

Specific  liquid  water  content  (kg  kg'1) 
Initial  NaCl  solution  molality 


Temperature 

CO 

0.1 

mol  kg~2 

0.01 
mol  kg"1 

0.001 
mol  kg'1 

-19.20 

0.0039 

-17.30 

0.0056 

-13.60 

0.0083 

-10.40 

0.0147 

0.0026 

-7.90 

0.0257 

0.0062 

-6.00 

0.0298 

0.0045 

-4.20 

0.0387 

0.0065 

-3.26 

0.0559 

0.0102 

-2.80 

0.0571 

0.0122 

-2.18 

0.0669 

0.0132 

-1.75 

0.0796 

0.0133 

-1.63 

0.0787 

0.0133 

-1.50 

0.0817 

0.0123 

-1.33 

0.0828 

0.0124 

-1.20 

0.0877 

0.0124 

-1.13 

0.0992 

0.0133 

-1.01 

0.1041 

0.0115 

-0.97 

0.1032 

0.0115 

-0.86 

0.1148 

0.0115 

-0.80 

0.1207 

0.0086 

-0.65 

0.1372 

0.0125 

-0.58 

0.1488 

0.0134 

-0.50 

0.1634 

0.0144 

-0.34 

0.1897 

0.0173 

-0.27 

0.2062 

0.0154 

-0.21 

0.0144 

-0.19 

0.0144 

0.0089 

-0.14 

0.0154 

0.0099 

-0.08 

0.0154 

0.0139 

-0.02 

0.0280 

0.0199 

Note:  The  equilibrating  solutions  of  the  pastes 
were  initially  0.1, 0.01,  and  0.001  mol  kg-1 
NaCl.  Data  are  not  presented  for  temperatures 
at  which  no  liquid  water  was  detected. 


resented  by  H20(cr,I)]  and  hydrohalite,  NaCl*2H20(cr)  (Archer  1992).  Our  calcu¬ 
lations  of  solution  molalities  assumed  implicitly  that  all  of  the  NaCl  remained  in 
solution  as  ice  was  formed.  We  did  not  calculate  the  molalities  of  solutions  for  obser¬ 
vations  below  252  K,  since  the  possible  formation  of  hydrohalite  at  these  tempera¬ 
tures  invalidated  our  assumption  that  all  of  the  NaCl  remained  in  the  liquid  phase. 
The  molality  of  the  liquid  solution  at  a  particular  temperature  of  interest  Tf  was  calcu¬ 
lated  by 


mTf  = 


Tf 


eT 


(22) 


where  wij  =  initial  NaCl  molality  of  solutions  of  pastes  (mol  kg-1) 

rtlff  =  NaCl  molality  of  solutions  of  pastes  at  temperature  Tf  (mol  kg-1) 
0j  =  initial  specific  liquid-water  content  (kg  kg-1) 

Qj(  =  specific  liquid-water  content  at  T  =  Tf  (kg  kg-1). 
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Specific  Mass  of  Unfrozen  Water  (kg  kg-1)  Specific  Mass  of  Unfrozen  Water  (kg  kg-1)  Specific  Mass  of  Unfrozen  Water  (kg  kg-i) 


Figure  1.  Unfrozen 
water  contents ,  as 
measured  by  pulsed 
NMR,  of  freezing  ka- 
olinite  pastes  cooled 
from  0X1  to  -66.6  X. 
The  equilibrating  so¬ 
lutions  of  the  pastes 
were  initially  0.1, 
0.01 ,  and  0.001  mol 
kg-1  NaCl. 


Figure  2.  Unfrozen 
water  contents,  as 
measured  by  pulsed 
NMR,  of  freezing  ka- 
olinite  pastes  warmed 
from  -66.6 X  to  OX. 
The  equilibrating  so¬ 
lutions  of  the  pastes 
were  initially  0.1, 
0.01,  and  0.001  mol 
kg-1  Nad 


Figure  3.  Unfrozen 
water  contents,  as 
measured  by  pulsed 
NMR,  of  freezing 
montmorillonite 
pastes  cooled  from 
OX  to  -66.6 X.  The 
equilibrating  solu¬ 
tions  of  the  pastes 
were  initially  0.1 , 
0.01,  and  0.001  mol 
kg-1  Nad 
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Specific  Mass  of  Unfrozen  Water  (kg  kg-i)  Specific  Mass  of  Unfrozen  Water  (kg  kg-i)  Specific  Mass  of  Unfrozen  Water  (kg  kg-i) 


Figure  4.  Unfrozen 
water  contents ,  as 
measured  by  pulsed 
NMR,  of  freezing 
montmorillonite 
pastes  warmed  from 
-66.6  V  to  OV.  The 
equilibrating  solu¬ 
tions  of  the  pastes 
were  initially  0.1 , 
0.01 ,  and  0.001  mol 
kg r1  Nad. 


Figure  5.  Unfrozen 
water  contents ,  as 
measured  by  pulsed 
NMR,  of  freezing 
sand  pastes  cooled 
from  OV  to -66.6  V. 
The  equilibrating  so¬ 
lutions  of  the  pastes 
were  initially  0.1 , 
0.01,  and  0.001  mol 
NaCl. 


Figure  6.  Unfrozen 
water  contents,  as 
measured  by  pulsed 
NMR,  of  freezing 
sand  pastes  warmed 
from  -66.6  VtoOV. 
The  equilibrating  so¬ 
lutions  of  the  pastes 
were  initially  0.1, 
0.01 ,  and  0.001  mol 
kg~J  NaCl . 


-60  -50  -40  -30  -20  -10  0 

Temperature  (°C) 
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The  solvent  and  solute  mole  fractions  could  be  calculated  by: 


(23) 

(24) 


*NaCl  " 


mTf 


mTf  + 


l 

Mh2o 


*H20  -  4  *NaCl 


where  Mh2o  is  the  molar  mass  of  water  (0.018  015  28  kg  mol  J). 


Freezing  point  of  the  solution  phase  in  bulk 

We  fitted  by  regression  the  freezing-point  depression  data  calculated  by  the  For¬ 
tran  program  FREZCHEM  (Marion  and  Grant  1994)  for  the  H20-NaCl  expression: 

T0  =  273 . 15  -  Kj mNaC1  (aq)  (25) 

where  Kf  is  an  approximate  cryoscopic  constant  for  the  H20-NaCl  system,  which 
we  found  to  be  4.12207  K  kg  mol-1.  This  compares  with  the  expected  cryoscopic 
constant  for  a  single  electrolyte  completely  disassociated  in  water  of  3.72  K  kg 
moH  (Atkins  1990). 


Molar  entropy  of  ice 

The  quantity  S^h2o  ®  Sm[H20(cr,I)]  at  a  temperature  and  pressure  of  interest 
can  be  calculated  by  adding  the  entropy  changes  due  to  cooling  and  freezing  to  the 
entropy  of  water  at  a  reference  point  (i.e.,  Tr  =  273.15  K  and  pr  =  0.1  MPa).  This 


quantity  is  calculated  by  evaluating 

Sm/Pr[H20(cr,I)]  =  S;/Tr<Pr[H20(l)] 

(26) 

+  ATTf7i-15S*m[H20(l),pr] 

(27) 

+  A1cr/IS;(H20,T  =  273.15,pr) 

(28) 

+  A^273'15S;[H20(cr,I),pr]. 

(29) 

The  starting  point  of  the  calculation  (eq  26)  is  the  molar  entropy  of  pure  water  at 
a  reference  temperature  ( Tr  =  273.15  K)  and  pressure  (pr  =  0.1  MPa),  for  which  the 
standard-state  molar  entropy,  S*mjr  /Pr  [H20(1)],  is  69.950  J  K_1  mol"1  (Chase  et  al. 
1985). 

Haida  et  al.  (1974)  reported  the  entropy  change  due  to  the  constant-pressure 
cooling  of  water  from  298.15  K  to  its  freezing  point,  At'-273  15S^1  [H20(1),  pr J  =6.615 
±  0.01  (eq  27).  This  quantity  can  also  be  calculated  with  the  equation-of-state  model 
for  water  (Hill  1990): 

%20(l)(rno/Prc)-SH20(l)(Trc/Prc)  =  -6-616  04  JK  1  mo1  *  (30) 

We  used  the  latter  value  in  our  calculations. 

The  entropy  change  associated  with  freezing  (eq  28),  A1cr  IS^(H20,  T  =  273.15, pr) 
is  calculated  by 


Acr,lSm(H20,r  =  273.15,pr)  = 


Acr,iHm(H2Q,T  =  273.15,pr) 
273.15 


(31) 


Haida  et  al.  (1974)  reported  a  molar  enthalpy  of  melting  of  6006.8  J  mol  1.  Giauque 
and  Stout  (1936)  earlier  reported  a  value  of  6007.0  ±  3.8  J  mol-1. 
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The  constant-pressure  entropy  change  due  to  changes  in  temperature  of  a  pure 
substance  (eq  29)  is  calculated  by 


An?  Sr 


T2  c 

=  p 


dr 


(32) 


As  with  many  solids,  the  heat  capacity  of  H20(cr,I)  as  a  function  of  temperature 
can  be  described  by  the  Maier-Kelly  (1932)  equation: 


Cp  =a+bT  +  ~2-  (33) 

Values  of  C*H2o(cr,i)  between  T  =  198.57  K  and  268.39  K  reported  by  Haida  et  al. 
(1974)  were  fitted  by  nonlinear  regression  (Marquardt  procedure)  to  eq  33  (SAS 
Institute  1985).  The  data  of  Haida  et  al.  (1974)  and  the  fitted  line  are  presented  in 
Figure  7.  The  regression  had  a  coefficient  of  determination  ( R 2)  of  0.999  997.  The 
following  parameter  estimates  (with  parameter-estimate  standard  errors  in  paren¬ 
theses)  were  obtained: 


a  =  -10.6644  (1.5999)  J  K_1  mol-1 
b  =  0.1698  (0.0046)  J  K-2  mol-1 
c  =  198  148.  (28  230.)  J  Kmol-1  . 
Equation  32  was  calculated  by 


Tj  Si dT  =  - HTno)]  +  KT(  -bTrw)+ 


(34) 


Molar  entropies  of  electrolyte  solutions 

The  procedure  for  calculating  the  molar  entropies  of  electrolyte  solutions  can  be 
derived  from  equations  found  in  standard  texts  on  chemical  thermodynamics  (e.g., 
p.  152-154  in  Lupis  [1983]).  The  entropy  of  a  mixture  or  solution  [SM  0  K-1)]  com¬ 
posed  of  nA  moles  of  A  (by  convention,  the  solvent)  and  nB  moles  of  B  (again,  by 
convention,  the  solute)  is 

5m  =5a  +  %  (35) 


Figure  7.  Constant- 
pressure  heat  capacity 
ofH20(cr,I)  under  p  = 
0.101325  MPa  andl  = 
200  K  to  268  K  as  re¬ 
ported  by  Haida  et  al. 
(1974).  The  line  was 
calculated  by  a  regres¬ 
sion  fit  of  the  data  to 
the  Maier-Kelly 
model. 
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Defining  the  total  number  of  moles  in  the  solution  as  nj  =  nA  +  nB,  this  is  equivalent 
to 


“  nA$m,A  +  nB^m,B 

where  Sm  A  is  the  molar  entropy  of  component  A  (J  K_1  moH)  and  Sm  M  is  the 
molar  entropy  of  the  mixture,  so 

=  =  -^A^ir^A  "*~*B®m,B*  (?7) 

The  calculation  of  Sm  M  is  facilitated  by  defining  the  molar  entropy  of  mixing, 
^mix^m,M'  as 

—  ^A^m,  A  ”  ^B^ir^B  (^) 

where  A  is  the  standard-state  molar  entropy  of  A.  The  molar  entropy  of  mixing 
for  the  mixture  is  also  equal  to  the  proportional  sum  of  the  molar  entropies  of  mix¬ 
ing  for  all  components  in  the  solution: 

^mix^m.M  “  *A^mix^m,A  +  *B^  mix^  m,B  (39) 

where  is  the  molar  entropy  of  mixing  for  component  A.  For  each  mixture 

component  i,  the  molar  entropy  of  mixing  can  be  calculated  by 


Arrux^nv  ~  ‘ 


^(AmixGm,;) 

=  -R  In  aX/ 

3T 

V 

'  l  3T  J 

(40) 


where  ax  z-  is  the  mole-fraction-based  activity  of  component  i.  So  that  the  molar  en¬ 
tropy  of  a  solution  can  be  calculated  by 


^m,M  —  —  XA  ^m,A  R^nax/A  RR 


^ain^A 
v  W 


+  xB 


Sm,B  “^ln ax,B  - RT\ - ^7 


d]nax 


(41) 


To  evaluate  this  equation,  the  following  quantities  must  be  determined  for  the  tem¬ 
peratures  of  interest: 

1.  The  standard-state  entropy  of  the  solvent,  water 

2.  The  standard-state  entropy  of  the  solute,  NaCl(aq) 

3.  The  activities  of  both  water  and  NaCl(aq) 

4.  The  temperature  derivatives  of  the  natural  logarithms  of  both  activities  of  the 
solvent  and  solute. 

Standard-state  entropy  of  liquid  water  at  subzero  temperatures.  The  standard-state  en¬ 
tropy  of  water  at  subzero  temperatures  is  estimated  with  the  heat  capacities  of 
supercooled  water. 

Speedy  (1987)  has  presented  equations  for  the  various  thermodynamic  proper- 
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Table  8.  Coefficients  to  Speedy's  (1987)  empirical  equations 
for  calculating  the  thermophysical  properties  of  supercooled 
water. 


X 

a  xlO3 
(K) 

Kj 

(bar-'1) 

(JK-imoH) 

cx 

-0.80 

20.0 

14.2 

Bx<°> 

1.802 180  3 

4.120 

25.952 

Bx(1) 

-0.941  698  0 

-1.130 

128.281 

Bx(2) 

0.905  507  0 

77.817 

14.2 

Bx(3) 

-0.80 

-78.143 

-221.405 

Bx(4) 

54.290 

-64.812 

max  resid 

1.2  ppm 

0.2% 

0.03% 

ties  of  supercooled  water.  His  expression  for  the  constant-pressure  heat  capacity  of 
supercooled  water  is 


4 

2 

n-0 


cp  =  5Xn)e"  +2C< 


c  e 


-1/2 


(42) 


where  and  Cq?  are  fitted  coefficients.  The  estimated  values  of  these 

coefficients  are  presented  in  Table  8.  The  parameter  e  is  a  reduced  temperature 


T-Ts' 

T 

1s 


,  Ts  being  a  limiting  temperature  that  is  assumed  to  be  exactly  227.15  K 


(^160C).  Figure  8  presents  the  constant-pressure  heat  capacity  of  water  as  calcu¬ 
lated  by  the  models  of  Hill  (1990)  and  Speedy  (1987). 

The  entropy  changes  of  supercooled  water  can  be  calculated  by  evaluating  the 
integral: 


Figure  8.  Constant-pressure  heat  capacity  ofH20(l)  and  H20(cr,I) 
under  p  =  0.101325  MPa  and  T  =  200K  to  300  K.  C v[H20(cr,I)] 
values  were  calculated  with  the  Maier-Kelly  model  CJH20(l)] 
values  between  T  =  200  K and  273  K  were  calculated  with  Speedy's 
(1987)  model  Cp[H2Od)]  values  between  T  =  274  K  and  300  K 
were  calculated  with  Hill's  (1990)  equation-of-state  model 
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T{  cp,h2o(1)  dT  Bc]M4  Tru'4) 

t  T  4  Tc4 


B^-4B^I7>3-7^3) 


3Tf 


+  ^ 


Bg»-3B®+6B«1(T,2-T„2) 


2t; 


Bg>  -2Bg>  4-3B0)  -4BW](r,  -T„) 


+  2Cc  arctan 


/  - \ 

5-i 

wTs  ; 


-  2Cr  arctan 


Lrw  _  | 

vj  Ts  ; 


[Bg)  -Bg>  +B® -BjP  +BW]x[ln(Tf)-ln(Tw)].  (43) 


Standard-state  entropy  ofNaCl(aq).  The  standard-state  entropy  of  NaCl(aq)  is  calcu¬ 
lated  similarly.  Initially,  we  can  define  the  apparent  molal  constant  pressure  heat 
capacity  of  the  solute,  ^Cp /NaC1(aq)  G  K-1  mol-1),  in  a  NaCl  aqueous  solution  as 


p,NaCl(aq) 


-  nH20(l)C®H20(l) 


MNaCl(aq) 


(44) 


where  C„  m  =  constant-pressure  molar  heat  capacity  of  the  solution 
Q  K-1  moF) 

C®h2o(1)  =  constant-pressure  heat  capacity  of  liquid  water  in  its 
standard  state  (i.e.,  pure)  (]  K-1  mol-1) 
nH20(i)  =  amount  of  water  in  the  solution  (mol) 
nNaCl(aq)  =  amount  of  NaCl  in  the  solution  (mol). 

The  change  in  entropy  for  the  solute  in  its  reference  state  is  calculated  from 


ASNaCl(aq)  =  l 


rn  6/^0 

1  Ip,  NaCl(aq) 


dr 


(45) 


The  apparent  molal  constant-pressure  heat  capacity  was  estimated  with  the  Pitzer 
model: 


,  C®H20(1)  _  Cp{mr)  ,  v|zNaza|Ac  ,_f  l  + 

■P,NaCl(aq)  =  £  MjWT  J 

“  2vNaVclRT2[(m-mr)BSaCl  +("j2  -wr)vNazNaCNaCl 


(46) 
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where  m  =  molal  concentration  of  the  solute  (mol  kg-1) 
nY  =  amount  of  NaCl  in  reference  solution  (mol) 
vNa  =  stoichiometric  number  of  Na  in  NaCl  (dimensionless) 

VC1  =  stoichiometric  number  of  Cl  in  NaCl  (dimensionless) 
v  =  vNa  +VC1  (dimensionless) 

Cp(mr)  =  constant-pressure  heat  capacity  of  a  solution  with  a  molality  equal 
to  mY  (J  K"1  mol-1) 

zNa  =  charge  number  of  the  Na+  cation  (dimensionless) 
zq  =  charge  number  of  the  Cl"  anion  (dimensionless) 

Ac  =  Debye-Hiickel  coefficient  for  apparent  molar  constant-pressure  heat 
capacity  (J  K_1  mol"1) 
b  =  constant  equal  to  1.2  (kg1/2  mol"1/2) 

Im  =  molality-based  ionic  strength  (mol  kg"1) 

Ir  =  molality-based  ionic  strength  of  the  reference  solution  (mol  kg-1). 

The  quantities  B^aC1  (kg  K"2  mol"1)  and  aci  (kS2  K_2  mol“2)  are  defined  by 


BNaC\  = 


NaCl 


dT2 


2 

+  — 
T 


ap(0) 


NaCl 


dT 


WaCll 

2 

+  — 

^P(NaCl 

dT2 

T 

dT 

L  J 

V 

V , 

and 


CNaCl  = 


rr-,2r(°)  1 

°  '-NaCl 

2 

+  — 

roc<01  1 

“'-NaCl 

+  < 

1 

°  '“NaCl 

2 

H - 

^NaCl 

dT2 

T 

dT 

dT2 

T 

dT 

L  J 

V 

- 

V 

L 

V 

L  J 

V . 

(47) 


(48) 


where  P^aC1  and  P^C1  are  ion-interaction  parameters  of  the  Pitzer  model  (kg 
mol-1);  C^cl  and  C^cl  are  ion-interaction  parameters  of  the  Pitzer  model  (kg2 
mol-2);  a  is  a  constant  in  the  Pitzer  model  (=  2.0  kg-1/2  mol-1/2);  and  a2  is  a  con¬ 
stant  in  the  Pitzer  model  (as  revised  by  Archer)  (=  2.0  kg-1/2  mol-1/2).  As  can  be 
inferred  from  eq47  and  48,  P($aC1,  P^aC1/  Cj^cl,  and  C^q  are  functions  of  T  and 
p.  Archer  (1992)  deveoped  lengthy  formulae  to  calculate  these  parameters. 

Activity  0/H2O.  The  activity  of  water,  uh20(1)'  111  an  aqueous  NaCl  solution  can  be 
calculated  by 

flH20(l)  =  exp[-<l)H20(l)vmNaCl(aq)]  (49) 

where  <t>H20(l)  is  the  osmotic  coefficient  of  water  (dimensionless)  (Stokes  1991). 
The  osmotic  coefficient  was  calculated  from  the  Pitzer  model  via 


18 


[fe2o(i)-i]=-|wciKY^^ 

+  m2VN^Va[p(Q)a+pWcl  expf-a^)] 

+  ffl2  4VNay  —  [c^cl  +  Cg>a  exp(-a2^) 


(50) 


Activity  ofNaCl(aq).  The  mole-fraction-based  activities  (see  App.  A)  that  are  called 
for  by  eq  41  are  related  to  the  molality-based  activities  typically  calculated  in  solu¬ 
tion  chemistry  by 

lnaxB  =  In  <1^3  +  ln(MAmBj  (51) 

where  MA  is  the  molar  mass  of  solvent  A  (kg  mol-1).  The  mean-ionic  activity  of 
NaCl(aq),  «+NaCl(aq)  (dimensionless)  is  calculated  by 


Y±NaCl(aq)m±NaCl(aq) 

fl±NaCl(aq)  = - © - 

™±NaCl(aq) 


(52) 


where  Y±NaCl(aq)  = 

m±NaCl(aq)  = 

m  ®  = 

±  NaCl(aq) 


mean-ionic  activity  coefficient  of  NaCl(aq)  (dimensionless) 
mean-ionic  molality  of  NaCl(aq)  (mol  kg-1) 

standard-state  mean-ionic  molality  of  NaCl(aq) 

(mol  kg-1). 


The  mean-ionic  activity  coefficient,  in  turn,  is  calculated  with  the  Pitzer  model 
for  a  mean-ionic  activity  coefficient  in  a  one-electrolyte  aqueous  solution: 


lnY±  HzNaZciK 


1  + 


+  m 


2vNavCl  2p(°)  |  ^NaCI 


l 


a2I„ 


1- 


2 

1  +  -  •^L  I  exp(-aTl^) 


+  m 


qr’(O)  +  4f'(1) 
^NaCl  +  ^NaCl 


2  4Vjvja  VC1 


(53) 


4r2  'N 


6  +  6a2^  +  3allm  +  3a%I%2 exp(-a2^) 


Temperature  derivatives  of  the  natural  logarithms  of  activities  of  both  H2OO)  and  NaCl(aq). 
Temperature  derivatives  were  calculated  numerically.  The  derivative  of  a  function 
/at  T  was  estimated  by  a  five-point  numerical  approximation, 

r{Tt)  =  ^[/(rf  - 25T) - 8/(T £  - 8T)  +  8/(Tf  +  S T)-f(T{  +  28T)]  (54) 

where  5 T  is  a  small  increment  chosen  to  be  0.1  K.  The  expected  error  of  this  esti¬ 
mate  should  be  on  the  order  of  5 T4  (Burden  and  Faires  1989). 
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Molar  volume  of  the  liquid  solution 

The  apparent  molar  volume  NaCl(aq),  *yNaci(aq)  (m3  niol-1),  in  a  simple  aque¬ 
ous  solution  is 


+V, 


0 

~  nu2o{i)vn2o(\) 

nNaCl(aq) 


NaCl(aq)  : 


(55) 


where  ^  VNaQ(aq)  =  apparent  molar  volume  of  NaCl(aq)  (m3  mol"1) 
yM  =  actual  solution  volume  (m3  mol-1) 

Vh20(i)  =  standard-state  molar  volume  of  water  (m3  mol"1). 

The  apparent  molar  volume  of  NaCl(aq)  was  calculated  by 

*  ^NaCl  =  ^NaCl  +  v|zNazCl  | +  b )  + 

2vNavCl^T[mNaClBNaCl  +  mNaCl(vNazCl)CNaCl]  (56) 

where  V^aC1  is  the  standard-state  molar  volume  of  NaCl  (m3  mol"1)  and  Ay  is  the 
Debye-Huckel  coefficient  for  apparent  molar  volume  (m3  kg"1/2  mol"1/2). 

The  quantities  B^aci  and  cNaCl  (both  kg  Pa-1  mol"1)  are  defined  by 


and 


rv  - 
LNaCl  - 


dCNaCl 

l 


\ 


( 


gCNaCl 

dp 


\ 


2|zNazCl|1/2 


(58) 


Temperature  derivative  of  the  solute  chemical  potential 

The  chemical  potential  of  the  solute  can  be  restated  in  terms  of  its  standard-state 
chemical  potential,  mean-ionic  molality,  and  mean-ionic  activity  coefficient: 


M'NaCl  =  PNaCl(aq) 

=  PNaCl(aq)  +  ^  T  In  fl±NaC1(aq) 

Y±NaCl(aq)m±NaCl(aq) 


=  M'NaCl(aq)  +  RT]n 


m±NaCl(aq) 


(59) 


where  j^NaCi(aq)  is  the  standard-state  chemical  potential  of  NaCl(aq)  (J  mol-1) 
and  fl±NaCi(aq)  1S  the  mean-ionic  activity  of  NaCl(aq)  (dimensionless).  The  tempera¬ 
ture  derivative  of  eq  59  was  calculated  numerically  by  eq  54. 
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Plots  of  capillary  pressures  against  liquid  specific  volumes 

The  plots  of  the  calculated  ice-solution  capillary  pressures  against  the  unfrozen- 
solution  specific  volumes  are  presented  in  the  following  figures.  In  general,  these 
plots  indicate  that  the  relationship  between  capillary  pressure  and  unfrozen  solu¬ 
tion  volume  was  unaffected  by  the  composition  of  the  liquid  phases. 

The  pc-  V  relationship  for  cooling  of  kaolinite  pastes  is  presented  in  Figure  9. 
The  general  trend  of  the  data  points  was  a  single  line  defining  the  pc  -  V  relation¬ 
ship  for  this  material.  For  the  kaolinite  pastes  that  were  initially  equilibrated  with 
0.1-mol  kg-1  NaCl  solutions,  however,  a  pronounced  deviation  from  this  line  was 
apparent  at  the  lowest  unfrozen  solution  contents.  Similar  deviations  were  found 
for  some  of  the  sand  and  montmorillonite  pastes,  but  only  for  the  pastes  washed 
with  0.1-mol  kg-1  NaCl  solutions.  This  is  apparently  a  systematic  error,  but  we 
have  been  unable  to  determine  its  source.  The  warming  curve  for  the  kaolinite 
pastes  is  presented  in  Figure  10.  No  hysteresis  was  apparent  when  comparing  the 
cooling  and  warming  curves  of  these  pastes,  though  here  again  there  is  an  appar¬ 
ent  systematic  error  in  some  of  the  data  collected  from  the  pastes  washed  with  0.1- 
mol  kg-1  NaCl  solutions. 

The  pc  -  V  relationships  for  cooling  and  warming  of  montmorillonite  pastes  are 
presented  in  Figures  11  and  12,  respectively.  Unlike  the  kaolinite  data,  the  general 
trend  traced  by  the  freezing  curves  for  montmorillonite  is  much  less  smoothly  curved 
than  for  kaolinite.  In  addition,  there  is  far  more  scatter  in  the  data.  We  believe  that 
the  abrupt  decrease  in  specific  volumes  at  higher  capillary  pressures  is  due  to  the 
large  relative  volume  of  liquid  that  was  held  in  the  interlamellar  spaces  of  this 
clay — in  these  narrow  spaces,  capillary  forces  play  a  less  significant  role  in  main¬ 
taining  the  water  unfrozen  than  does  the  so-called  'surface  melting  of  the  ice'  (Brun 
et  al.  1977,  Dash  1989). 

The  pc-  V relationships  for  cooling  and  warming  of  sand  pastes  are  presented  in 
Figures  13  and  14,  respectively.  As  with  montmorillonite,  the  general  trend  traced 
by  the  freezing  curves  of  this  clay  is  much  less  smoothly  curved  than  for  kaolinite. 
We  believe  that  this  abrupt  transition  is  due  to  the  coarse  and  fairly  homogeneous 
pore-size  distribution  of  these  pastes,  which  allowed  the  solutions  to  freeze  at  higher 
and  generally  more  uniform  temperature. 


Figure  9.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  kaolinite  pastes  cooled  from  0X1  to 
- 66.6  X.  The  equilibrating  solutions  of  the  pastes  were  initially  0.1 , 0.01, 
and  0.001  mol  kg"1  NaCl. 
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Figure  10.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  kaolinite  pastes  warmed  from  -66.6  X 
to  0X1.  The  equilibrating  solutions  of  the  pastes  were  initially  0.1 ,  0.01, 
and  0.001  mol  kg r_1  NaCL 


Figure  11.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  montmorillonite  pastes  cooled  from 
OX  to  -66.6  X.  The  equilibrating  solutions  of  the  pastes  were  initially 
0.1,  0.01,  and  0.001  mol  kg-1  NaCL 
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Figure  12.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  montmorillonite  pastes  warmed  from 
-66.6X1  to  OX!.  The  equilibrating  solutions  of  the  pastes  were  initially 
0.1,  0.01 ,  and  0.001  mol  kg'1  NaCl. 


Figure  13.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  sand  pastes  cooled  from  0X2  to 
-66.6  ‘C.  The  equilibrating  solutions  of  the  pastes  were  initially  0.1 , 0.01, 
and  0.001  mol  kjf1  NaCl 
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Figure  14.  Relationships  between  unfrozen-solution  specific  volumes  and 
ice-solution  capillary  pressures  for  sand  pastes  warmed  from  -66.6  T2  to 
0  U.  The  equilibrating  solutions  of  the  pastes  were  initially  0.1, 0.01 ,  and 
0.001  mol  kg ^  NaCl. 


CONCLUSION 

The  minerals  analyzed  in  this  study  may  be  thought  to  represent  the  range  of  ap¬ 
plicability  of  capillary  theory  to  explain  the  freezing  behavior  of  porous  solids.  The 
kaolinite  samples  were  not  composed  of  an  expanding-lattice  clay  with  intermedi¬ 
ate  particle  size,  and  the  experimental  results  from  this  material  supported  the  theory 
well.  Sand  and  montmorillonite  did  not  support  it  as  well:  sand  because  its  particle 
size  (and,  therefore,  pore  radii)  was  too  large,  and  montmorillonite  because  it  has  a 
much  smaller  particle  size  and  because  relatively  large  volumes  of  water  can  be 
held  in  its  interlamellar  spaces  for  which  capillary  forces  are  less  important.  Our 
study  did  indicate  that  capillarity  theory  coupled  with  the  Pitzer  model  describes 
acceptably  the  effect  of  an  electrolyte  on  the  freezing  curves  of  minerals  of  interme¬ 
diate  size  between  sand  and  montmorillonite. 
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APPENDIX  A:  EFFECT  OF  CONCENTRATION  SCALES  ON  ACTIVITIES 


The  value  of  an  activity  of  a  solute  depends  on  the  concentration  scale  that  is 
chosen.  This  is  because  the  standard-state  chemical  potentials  for  different  refer¬ 
ence  functions  are  themselves  different  (Mills  1988).  For  a  molality  reference  func¬ 
tion,  the  standard-state  chemical  potential  for  solute  B  is 


M-m,  b(t)  ~ 


f 

p(T,pe,mB  )-RT In  -^g- 
V  ’  \mB . 


(60) 


and  the  standard-state  chemical  potential  for  solute  B  for  a  mole-fraction  reference 
function  is 


(T)  =  [|t(T ,  pG ,  nrB,. . . )  -  RT  ln(xB  )j°° . 
Subtracting  eq  60  from  eq  61  yields 


This  is  equivalent  to 


(T)-P®B(T)  =  -RT 


ln(xB)-ln 


VmB  j 


In 


f  \ 


1 


lMA 


+  mB 


-In 


my 

mf 


(61) 


(62) 


(63) 


and 


H?,b(T)-VL%b(T)  =  -RT 


In 


0  , 
mB  J 


(  \ 

mB 

-In 

mB 

TJ-+WB 

Ma 

lmB  J 

which  results  in 

H?b(T)  =  bCO  -  RTln(MAmf ). 
The  molality-based  activity  of  solute  B  is 


V  B  =  exp 


P-B 

RT 


and  its  mole-fraction-based  activity  is 


ax,B  =  exP 


hB  ~  P®,B 


RT 


(64) 


(65) 


(66) 


(67) 
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can  derive  immediately 


lnflX/B  =  lnVB  +  lnl 
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